Dynamical Casimir Effect in two-atom cavity QED 
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We study analytically and numerically the dynamical Casimir effect in a cavity containing two 
stationary 2-level atoms that interact with the resonance field mode via the Tavis-Cummings Hamil- 
tonian. We determine the modulation frequencies for which the field and atomic excitations are 
generated and study the corresponding dynamical behaviors in the absence of damping. It is shown 
that the two-atom setup allows for monitoring of photon generation without interrupting the growth, 
and different entangled states can be generated during the process. 



(N 

o 

00 



I 



> 

cn 

rn 
O 
(N 



% 



PACS numbers: 42.50.Pq, 32.80.-t, 42.50.Ct, 42.50.Hz 

Introduction.-ln view of the recent progress [Ij in ex- 
periments on the observation of the so called Dynamical 
Casimir Effect (DCE) [2| , the problem of detecting pho- 
tons generated from the initial vacuum state becomes 
quite actual. It was shown long ago Q that the presence 
of a detector can change significantly the statistics (in- 
cluding the mean number) of created quanta, compared 
with the predictions made for an idealized empty cav- 
ity model. Therefore it is necessary to study in detail 
different detection schemes. At least two main schemes 
were proposed until now. In the so called MIR experi- 
ment the quanta of the microwave field are to be detected 
by an antenna put inside the closed cavity [J|. Another 
idea was to use as detectors real Rydberg atoms gass- 
ing through the cavity [3, 0-0| or "artificial atoms" [8J in 
the case of Circuit QED systems, such as those described 
in [il, 0| • The simplest solutions for a detector modeled 
as a single two- level atom were obtained in [3J, |5|, and 
recently results of more detailed theoretical and numeri- 
cal studies of the atom-field interaction during the DCE 
were presented in iQnljJl . Three- level models of detectors 
were considered in [i3,ll3- It seems important to analyze 
different configurations to choose the optimal scheme. 

Here we study how the DCE dynamics is affected by 
the presence of two 2-level atoms (detectors) interacting 
with a single resonance cavity field mode. Our starting 
point is the Hamiltonian (we set h = I) 



Ha = uJtn + Y^ [-^cr| +.9j(acr+ +aV^. ) -ixt{a 



,t2 



where a (a^) is the cavity annihilation (creation) oper- 
ator and n = a'^a is the photon number operator. The 
Pauli operators are defined as (t| = |ej)(ej| — \gj){gj\, 
^J = lffi)(ej|, (7J = |ej)(5j|, where |g^) and \ej) are the 
ground and excited states of the j-th atom (j = 1,2), 
respectively, ilj and gj are the atomic transition fre- 
quencies and the atom-field coupling constants (assumed 
real for simplicity). If Xt = Oj then Hq is the spe- 
cial case of the known Tavis-Cummings Hamiltonian [14l | 
studied in numerous papers (see, e.g., |15l - [l8l | and refer- 
ences therein) . Physical realizations of this Hamiltonian 
(which holds for \gj\ <^ D,j) were demonstrated in [1^ for 
trapped ions and in [201 for the Circuit QED systems. 
The last term in Hq describes the effect of photon cre- 



ation (equivalent to squeezing) in a cavity whose fun- 
damental eigenfrequency varies in time due to the mo- 
tion of a boundary [s], l2l|, |2^. We suppose that the 
boundary performs harmonic oscillations at the modu- 
lation frequency rj. Then the instantaneous cavity eigen- 
frequency depends on time as Wj = wq +£ sin{r]t), where e 
is the small modulation amplitude. Normalizing the un- 
perturbed cavity frequency to wg = 1, we write the modu- 
lation frequency as 77 = 2 (1 -I- x), where a; is a small reso- 
nance shift. For a weak modulation, \e\ ^ 1, we can write 
to the first order in e: Xt = {4uJt)^^dut/dt ~ 2qcos{rit) 
[3, I21I, [22], where q = e {1 + x) /4. Moreover, the term 
LOtn in Hq can be replaced simply by n, as soon as the 
main effect of modulation is due to the presence of op- 
erators a^ and a^'^ in the squeezing part of Hq, but not 
due to the photon number preserving part tuta^a. 

In the empty cavity, the resonance generation of many 
photons is achieved for x = (being impossible if \x\ > 
\e\ 23]). On the other hand, it was shown [3] that no 
more than two photons can be created in the presence 
of a single atom if \e\ <C |5i|, and this can happen if 
\x\ ~ \gi\. Our aim is to find the resonance regimes in 
the presence of two atoms for different relations between 
the parameters e, gj and Qj. We show that there are 
two types of resonances. For some distinguished values 
of a; 7^ at most two photons can be created. But under 
certain conditions, the multiphoton generation becomes 
possible again for x ~ (contrary to the one-atom case), 
even if |e| ^ \gi\. This interesting result is one of the 
main motivations for this publication. 

The dynamics of the closed system (atoms -I- field 
mode) is governed (neglecting dissipation) by the 
Schrodinger equation id\'^{t))/dt = iJo|*(i))- To find 
analytical solutions we go to the interaction picture: 
|*(t)) = exp[-it{ri/2){n + af/2 + al/2)]\'ilj{t)), since 
the Hamiltonian acting upon the new wavefunction \ip{t)) 
becomes time independent after the Rotating Wave Ap- 
proximation (RWA): 






j=i 



-2-"^- 



iq{a —a' )—xn, 



where A^ = 1 — Qj . We expand the wavefunction in the 
atom and Fock bases as follows: 



m)) = E 



m=0 



«™We-*(2-+^^+^^)*/2|ffi)|52)k)+6„We-'(^^-^^)*/2|5i)|e2)|m) 



Then the Sehrodinger equation with Hamihonian Hj leads to the set of coupled differential equations 



Cm-1 = -i5i\/ma„e"'^'* - ig2^m ~ ldm-2e^^^* + qWm-iCm-i 

dm-2 = "igWm - 16™_ie"*^i* - ig2^m - lc„j_ie"*^^* + (7W„-2rfm-2, 



'igi\/m~- ld„i_2e*^'* - ig2\/ma„ie ^'^^'^ + qWm-ibm-i 



(1) 



(2) 
(3) 
(4) 
(5) 



where WmO„ 



Vm(m-l)0„_2e-2«* _ ^(m + 1) (m + 2)0,„+2e'"*. 



Weak modulation with atoms in resonance. -This 
regime is defined by the inequality |e| <C G = i/ffi + g| ■ 
If two atoms are in resonance, Ai = A2 = 0, the solu- 
tion to Eqs. (HI-® in the absence of external modulation 
{q = 0) is (for m > 2) 



Q,/9 = + : 



J-;^^exp(mGLli), 



J'm-2 — 



E 

a.p=+,- 

G 



Kf,J-;^^exp(azGLf„i), 



(6) 

(7) 



Q,^ = 



where J^^'^ are constant coefficients, 
1 =F 2i?m _ 2.9132 



E ai^,^;?,^exp(mGL^i) 



"1] 



(8) 
(9) 



V^ = 



Rn 



Ip^Jm (ra — 1) G 

iyiTv^^r^^^r^, L± = V?" - 1/2 ± i?« 



Substituting now expressions ([B])-® back into Eqs. ([2])- 
([5]) and assuming that T^ are slowly varying func- 
tions of time, one can verify that for specific values of 
the resonance shift x some of these functions become 
multiplied by imaginary exponentials with large argu- 
ments (compared to g), while others are multiplied by 
time-independent coefficients, so one is allowed to per- 
form the RWA and obtain simplified effective dynam- 
ics. We find that for the initial zero-excitation state 
I5i)l52)|0) 9.t most two photons can be created when- 
ever G 1^4 — L2 I ^ q- The resonant regimes occur for 
1x — —aGL2 (with a, /3 = +, — ), when the only nonzero 
amplitudes (neglecting small terms of the order of e/G) 
are Oq = cos{qtRf-j) (it does not depend on the sign of a) 



and ^"2" = Rp sm{qtRfi)/V^, where R± = i J2 ± R. 



For a single atom (32 
that R+ = 1 and -R_ = 



0) one has Rm = 1/2, so 
Then the only resonances 



with a periodic creation of at most two photons happen 
for X = ±\gi\/V2. In this case oq = cos{qt), while the 
only other nonzero coefficients are 7^2^^ — sm{qt)/\/2 in 
accordance with Q . In the presence of the second atom, 
new resonances become possible. If \g2\ <C \gi\, then 
these additional resonance frequencies have the values 
X w ±|5i|/2. However, since i?_ w p\/2 ^ 1 in this 
case, the corresponding dynamics is quite slow and the 
probability of the photon creation is small, too. 

The most interesting situation takes place if \gi\ = 
\g2\- Then Rm = m — 1/2 and L~ = 0. We still have 
the resonances at a; = ±|3i|\/3/2, when no more than 
two photons can be created from the initial ground state, 
since the only nonzero coefficients in this case are oq — 
cos(\/2/3gi) and 7-2^^ — r sin(-y/2/3qi)/-\/3, where r = 
92/ gi — ±1- But two other resonances merge in the 
single one at x = 0. In this case, solving Eqs. (HI)-® 
with g = 0, one can write (for to > 2) 

am = r [WraEmit] + XmE+{t) + y„,] , 

6„_i = Vl - (2"i)-i [W™S™(t) - XmE+{t)] + Zm, 

2to 



Mn—2 ^^n 



TO, — 1 



1 



m y^m{m- 1) 



yr, 



where E^it) — exp[±i5i-y/2(2m — l)t]. In the presence 
of additional terms proportional to the small parameter 
q <C G in Eqs. dH-®, the coefficients W™, X,n, ym and 
Zm become time-dependent. For the standard atomless 
DCE resonance ?/ = 2, assuming that |>Vm| , \Xm\ ^ 1 
for all TO, we perform the RWA and find that Zm{t) — 0, 
meaning that 6m(i),c„i(t) « for all times. Only func- 
tions y^n vary slowly with time according to the equations 



yn 



/— —2m -3 

y/m{m - 1) -ym-2 

Zm — 1 



r, TT r TO, — 1 2to + 1 

- V(™+l)(™ + 2)^— -yrn+2 

771-1-1 Zm — 1 



(10) 



with the initial condition 3^m(0) — rSmo- Therefore even- 
tually all (even) coefficients 3^™ become different from 
zero, so that many photons can be created from the initial 
vacuum state. Eq. (|10p has two remarkable properties. 
First, it does not contain the atomic coupling coefficients. 
Second, the fractions in its right-hand side tend to the 
unit values for ?7i 3> 1, and in this limit Eq. (|10l) has the 
same form as the equation governing the evolution of the 
field amplitudes (in the Fock basis) in the empty cavity. 
Since the main contribution to the mean photon num- 
ber (n(t)) is given by the coefficients 3^^ with m >• 1 if 
{n) ^ 1, we can expect that after some transient time the 
photons will be steadily created with the same asymp- 
totical rate d\n{{n)) / d{et) as in the empty cavity. More- 
over, since |dm_2(i)| = [m/{m—l)] \am{t)\ , both atoms 
become excited simultaneously. Numerical calculations 
confirm these predictions, as shown in Fig. [U where we 
plot the mean photon number (n) and the probability of 
double excitation P(ei,e2} for parameters gi — A x 10^^ 
and e = 2 X 10^"^ |26|]. Part (a) shows the role of the de- 
tuning parameter x when 52 = 9i'- the photon creation 
and atomic excitations practically stop for x ^ e. Part 
(b) shows the influence of disbalance 32 — 5i when a; = 0: 
again, all effects practically disappear if |g2 ^ gi\ ^ s. 

The mean number of photons for x = is smaller than 
that in the empty-cavity case, (f^o(i)) ~ sinh (ei/2), due 
to initial transient processes, when the atomic popula- 
tions attain stationary values: one can see that the line 
{n{t)) can be obtained from {no(t)) by some positive shift 
in time. Therefore the a; = resonance for \gi\ = \g2\ is 
interesting from the point of view of detecting Casimir 
photons, since the atoms get excited simultaneously with- 
out interrupting the photon generation process. 

If the second atom is in the dispersive regime, I52I "^ 
IA2I (while Ai = 0), we define the dispersive shift 
62 = g\l ^2 and repeating the previous steps we find that 
for \82\ ^ I51I the photon generation occurs for the reso- 
nance shifts 2x = (3/2) 82 ± G2 with G2 = x/2gl + Sj/i. 
The resulting nonzero probability amplitudes read: ao = 

cos (^qt^/T±h/{2Ch 



02 



Cl 



_ ~i{3/2)52t 



[We 



-lG2t 



Xe 



iG2t-\ 



{W[l-52/i2G2)]e 



-iG2t 



V2gi 

-A'[l + <52/(2G2)]e^^^*}, 



61 ~ V2{g2/A2)e' 



■iA2t 



02, 



(g2/A2)e-'^^'ci, 



y/l±S2/{2G2) . ( ^ . ,^ ,,,„, 

-= sm ( qtV 1 ± 02/(262) 



At most two photons can be created in this case. 

Dispersive regimes-Many photons can be generated 
from vacuum if both atoms are in the dispersive regime, 
\gj\ <C |Aj|. In this case, instead of solving coupled dif- 
ferential equations it is convenient to write the wavefunc- 
tion IV'(i)) as [ll|,[ll |V'(i)) = C/'f exp (-ii7e/i) C/|^(0)) , 
where the effective Hamiltonian Hef = UHjU^ is defined 




FIG. 1: (Color online) The mean photon number (solid lines) 
and atomic excitation probabilities (dashed lines) as functions 
of dimensionless time et. (a) The influence of nonzero detun- 
ing X for gi = g2- (b) The influence of disbalance g2 — gi for 
X — 0. Numerical values of parameters are given in the text. 



by means of the unitary operator U — exp(y). Choos- 
ing y = a^ (C20'2^ + Cio'r) ^ '^■C- (where Q — gj/^j are 
small parameters, |Cj| <C 1, ^ = 1,2) and expanding the 
exponentials in Taylor's series we obtain to the second 
order in (j [assuming ©(Ci) ^ ^'(^2)] 



Hef - 



2 ^. 
-(x + Sial + 52al)n - X! ~^ 



3 ^« 



J = l 



-C1C2 



^i "■2 ~ 2.iqaja2 + /i.e. 



iq[{l + (:lal+Qlal)a'-h.c] 



(11) 



Here 5j = 9]l^3 ^'^'^ ll^*^ dispersive shifts (j — 1,2). In 
view of the perturbative expansion the effective Hamil- 
tonian (|TT1l is valid roughly for times |i5i|i ^ 1. For 
(7 = it describes the indirect interaction between the two 



atoms via the cavity field [2J|. Since the state |(?i)|52)|0) 
is the common eigenstate of Y and is'\'^2 + /i-c) with 
nufl eigenvalues, one has t/|5i)|52)|0) = I5i)l52)|0>: so 
the term (crj'crj -l- h.c^ can be dropped out if |V-'(0)) = 
I5i)l52)|0)- Besides, if the coefficient in front of n in 
Eq. (TTl]) is adjusted to zero, the photon generation 
term iq{\ + Ci ""i + Cl""!) '^^^ becomes resonant, while 
the term 2iqC^\C^2{p\ (^2 ~^-^-) ^^"^ be neglected for initial 
times. In this case the wavefunction for the resonance 
shift X = 8\ + 82 reads as (neglecting a global phase) 
IV'(i)) = f7l'A(l - C^) 151)1^2) |0), where the squeezing 
operator A(u) = &^\i\u qtia^'^ — a^)\ has the property 
[lllll k^{v)ak{y) =C^a + S^a^, with C„ = cosh(2t;gi), 
Sy = sinh {2vqt), and C^ 



c? + cl- 



Average values of the main observable quantities are 
as follows (to the second order in Q): 

{n{t)) = (1 - C') sinh^ [2gi (1 - C')] , 
Peiit) - Cf{n{t)), P,2{t) ^ Ci {n{t)) , 

\2\ 



1 



((A^±n = ^{C^ + (1-C^)exp[±4,i(l-C^)]} 



where X+ 



/V2 and X_ 



/{V2t 



are the field quadratures. Moreover, for times |(5i| i <C 1 
the probability P{ei,e2} of detecting simultaneously both 
atoms in their excited states is proportional to Ci j so it 
is very small. Therefore, by measuring Pel or Pe2 one 
can estimate the mean photon number. In Fig. [5^ we 
show the behavior of (n), P^i, Pe2 and P{ei,e2} for pa- 
rameters e = 2 X 10"^, 51 = 4 X 10"^, 32 = 3 x 10"^, 
Ai — lOgi, A2 — 15^2, and x = di + 62- We see that 
many photons are created and the atomic populations 
are proportional to the mean photon number, while the 
probability of double atomic excitation is very small. 



If 



EL(A,+35,) 



3> q and the resonance shift is 



tuned to 2a: = — 22j=i i^j + ^j) ^'^^^ ^1 ^ — A2, then 
the photon generation term becomes off-resonant and the 
only resonant term 2iq(!^i(2{<^i<^2 ~ h.c.) survives in the 
interaction part of the effective Hamiltonian (|lip even to 
higher orders in Ci , extending its validity beyond the pre- 
vious condition |(5i|t <C 1. In this case only the atomic 
excitations are generated at a rather small rate 2qC,iC,2 
and the probability of detecting both atoms simultane- 
ously in the excited states is (l — Ci ~ Cl) sin^(2(7tCi(^2)- 
In Fig. ^Bp we show the behavior of (n) and P[ei^e2} for 
parameters 51 = 4 x 10"^, 52 = 3 x 10"^, Ai = 0.22, 



A2 = -0.2, e = 2 X 10-3 and 2x = -E^=i (Aj +(5j), 
where we see that double atomic excitations are created 
while essentially the field remains in the vacuum state. 




FIG. 2: (Color online) The mean photon number and 
atomic excitation probabilities versus the dimensionless time 
et in the dispersive regimes: (a) x = 5i -\- 82', (b) 2x = 
— X]j=i (^i + ^j)- Other parameters are specified in the text. 



Other regimes.-li atom 1 is resonant (Ai = 0) and 
weakly coupled to the field {\gi\ <C s), while atom 2 is 
in the dispersive regime (I52I ^ 1^2 1), then we make 
the transformation with Y — a'' (C2'^2 ~^ '^^I'^t) ~ h.c, 
Ci — .91/(29) and C2 — <?2/A2. For the resonance shift 
X — 62 the effective Hamiltonian describing parametric 
amplification reads (after RWA) 



H. 



26n 



e/ 



^2 - k [(1 + ^I'^l + Cl<^2 



h.c] 



-<52 (1 - 2e?af + ^2^ 



n - 



Ic 



2^? 



fTi 



For the initial state |ei)|g2)|0) one has f/|ei)|g2)|0) = 
\e^)\g,)\0),som))=WA{l+ei-a)\ei)\g2m{npto 



a global phase). This yields the following average values: 



((AX 



n{t)) 
Pal 

2\ 



±) 



(l-ei-C2)sinh'[29(l+e?-C2)^], 

ei{n{t)), Fe2=C2< 



n{t)) , 

„±4g(l+4?^ 



ci)t 



where Pgi is the ground state probability of atom 1. Be- 
sides, the probability P{gi^e2} of finding simultaneously 
atom 1 in the ground state and atom 2 in the excited 
state is zero [to the second order in 0{£,i), 0{(2)]- 

Analogously, if both atoms are weakly coupled to the 
field, G <C l^l, then by performing the transformation 
with Y = ia^ {S.i^'^ + C2f^) — h.c. and £,j = gj/{2q) one 
obtains for x = Ai = A2 = the effective Hamiltonian 
[to the second order in £^j, for 0(^1) ^ 0(^2)] 



H,f=tq[{l 



e2„2 

^lo^l 



- ^2^2^ 



42 



- 2^i£_2(7ta+ - h.c] 



In these cases many photons can be created as well, and 
the atoms may serve to monitor the photon generation. 
Conclusions. -^e found that the two-atom nonstation- 
ary cavity QED is attractive from the point of view of 
producing different types of entangled states and detect- 
ing the DCE, because in specific regimes the atoms can 
acquire independent information about the field state 
without inhibiting the photon generation process. In par- 
ticular, we showed that in the realistic case when the ex- 
ternal modulation amplitude is much smaller than the 
atom-cavity coupling strengths, many photons, as well 
as atomic excitations, can be generated from the initial 
zero-excitation state even if both atoms are resonant with 
the unperturbed cavity field, contrary to the single 2- 
level atom scenario. Moreover, simply by adjusting the 
modulation frequency, keeping the other parameters un- 
altered, one can achieve the regime in which at most two 
photons are generated. If the atoms are off-resonant, 
then for the zero-excitation initial state many photons 
can be created for a specific modulation frequency; yet 
by appropriately tuning the modulation frequency one 
can achieve the regime in which only atomic excitations 
are generated. Furthermore, one can explore the regime 
in which one atom is resonant but weakly coupled to the 
field, while the other atom is in the dispersive regime - 



in this case many photons can be created from vacuum 
and the atoms monitor independently the process. This 
variety of possibihties can be useful for choosing optimal 
schemes of detecting the Casimir photons. In view of the 
results obtained, generalizations to the systems of three 
and more atoms could be quite interesting. But we leave 
this problem for another study. 
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